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ABSTRACT 

Context. Dark gamma-ray bursts (GRBs) are sources with a low optical-to-X-ray flux ratio. Proposed explanations for this darkness 
are: i) the GRB is at high redshift ii) dust in the GRB host galaxy absorbs the optical/NIR flux iii) GRBs have an intrinsically faint 
afterglow emission. 

Aims. We study two dark GRBs discovered by Swift, namely, GRB 100614A and GRB 100615A. These sources are bright in the 
X-rays, but no optical/NIR afterglow has been detected for either source, despite the efforts of several follow-up campaigns that have 
been performed since soon after the GRB explosion. 

Methods. We analyze the X-ray data and collect all the optical/NIR upper limits in literature for these bursts. We then build optical-to- 
X-ray spectral energy distributions (SEDs) at the times at which the reddest upper limits are available, and we model our SEDs with 
the attenuation curves of the Milky Way (MW), Small Magellanic Cloud (SMC), and one obtained for a sample of starburst galaxies. 
Results. We find that to explain the deepest NIR upper limits assuming either a MW or SMC extinction law, the visual extinction 
towards GRB 100614A is A v > 47 mag, while for GRB 100615A we obtain A v > 58 mag using data taken within one day after the 
burst and A v > 22 mag even 9.2 days after the trigger. 

Conclusions. A possible explanation to these unlikely values is that optical radiation and X-rays are not part of the same synchrotron 
spectrum. An alternative, or complementary explanation of the previous possibility, involves greyer extinction laws. A starburst 
attenuation curve gives Ay > 11 (A v > 15) for GRB 100614A (GRB 100615A) before 1 day after the burst, which is less extreme, 
despite still very high. Assuming high redshift in addition to extinction, implies that A v > 10 at z = 2 and A v > 4 - 5 at z = 5, 
regardless of the adopted extinction recipe. These lower limits are well above the A v computed for previous GRBs at known redshift, 
but not unlikely. A different, exotic possibility would be an extremely high redshift origin (z > 17 given the missing K detections). 
PopIII stars are expected to emerge at z ~ 20 and can produce GRBs with energies well above those inferred for our GRBs at these 
redshifts. However, high Nh values (above the Galactic ones) towards our GRBs challenge this scenario. Mid- and far-IR late afterglow 
(> 10 5 s after trigger) observations of these extreme class of GRBs can help us to differentiate between the proposed scenarios. 
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1. Introduction 

Long-duration gamma-ray bursts (GRBs) are high energy phe- 
nomena linked to the death of massive stars, emitting most of 
their radiation in the hundreds of keV range. The gamma-ray (or 
prompt) event is followed by an afterglow at longer wavelengths, 
which is crucial to understand the physics of these sources and 
investigate the nature of their surrounding medium. 

While the X-ray (0.1-10 keV) afterglow is virtually detected 
for all GRBs, the optical/near infrared (NIR) one is more elu- 
sive. The first optical afterglow detection (van Paradijs et al. 
1997) suggested the idea that every GRB had an optical coun- 
terpart. However, just a few months later, the search for the op- 
tical afterglow associated with GRB 970828 was unsuccessful 
(Groot et al. 1998), leading to the definition of 'dark burst' to 
underline a GRB with an X-ray counterpart, but not an optical 
one. GRB 970828 was not an isolated event, and the low suc- 
cess rate in the optical/NIR detection of afterglows became a hot 
GRB topic (see e.g., Fynbo et al. 2001 and Lazzati et al. 2002). 
Initially, this lack of detection was widely associated with the 
delay time between the GRB and the optical observations, since 
ground-based facilities could be on target only several hours af- 
ter the trigger, when the afterglow faded below their sensitivity 
threshold. 



The main scientific driver of the Swift satellite (Gehrels et al. 
2004) was to facilitate the GRB afterglow detection through a 
quick repointing with its narrow field instruments (XRT in the 
0.3-10 keV band, Burrows et al. 2005 and UVOT at UV/optical 
wavelengths, Roming et al. 2005) and a fast dissemination of 
the GRB coordinates worldwide. Despite these instruments re- 
point the target fewer than two minutes after the prompt event, 
the UVOT detection rate is just ~ 40% (Roming & Mason 
2006). However, the quick detection with XRT and the increas- 
ing number of ground-based automated facilities dramatically 
improved the optical follow-up success, reducing the fraction of 
dark GRBs. 

The optical darkness can be ascribed to different factors (see 
e.g., Fynbo et al. 2001, Perley et al. 2009). First, the GRB can be 
at high redshift, so that the Lyman a absorption prevents optical 
identifications. Second, dust in the GRB host galaxy or along 
the line of sight can absorb the optical afterglow counterpart. 
Finally, the optical faintness can have an intrinsic origin. 

The definition of dark GRB evolved as more data became 
available. Originally, a dark GRB was an event with an X-ray 
afterglow but no optical detection (Fynbo et al. 2001). Then, 
brightness and time limits were added to make this definition 
more specific, e.g., R > 23 within 12 hr after the prompt event. 
Finally, the physics of the GRB was involved in the dark/bright 
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dichotomy. The basic prediction of the fireball model (Meszaros 
& Rees 1997), which is commonly invoked to explain the after- 
glow, is that synchrotron radiation is responsible for the optical- 
to-X-ray emission. According to this model, the spectral index in 
the optical and X-rays is a function of the power-law index p of 
the electron energy distribution, and is (p - l)/2 or p/2 depend- 
ing where the synchrotron cooling (injection) frequency in the 
slow (fast) cooling regime is located with respect to the observed 
band (Sari et al. 1998). The simplest blast wave models (Wijers 
et al. 1997; Sari et al. 1998) predict p > 2, which translates to a 
spectral index > 0.5. In this framework, Jakobsson et al. (2004) 
proposed a classification in which dark GRBs had an optical-to- 
X-ray spectral index fiox < 0.5. Despite many GRBs behaving 
this way, it has been shown that p < 2 is achievable both obser- 
vationally and theoretically, by introducing a high-energy cut-off 
in the electron distribution (see e.g., Dai & Cheng 2001, Starling 
et al. 2008 and references therein). Thus, fiox can in principle 
be lower than 0.5 in the synchrotron model. Rol et al. (2005) 
proposed a classification criterion based on the faintest allowed 
extrapolation of the X-ray flux to the optical regime, inferring 
the p value from both the X-ray spectral and temporal indices 
according to the fireball closure relationships (see e.g. Sari et al. 
1998). 

Another approach was proposed by van der Horst et al. 
(2009). They noted that, regardless of many assumptions about 
the specific electron energy distribution, if both the optical and 
X-ray radiation are produced by synchrotron emission from the 
same source, the spectral indices in the optical (J3o) and X- 
ray (fix) bands are linked. In particular, fio = fix - 0.5 if the 
cooling frequency lies between the optical and the X-rays, and 
fio = fix otherwise. Thus, the optical-to-X-ray spectral index al- 
lowed range is fix - 0.5 < fiox ^ fix, with fiox - fix - 0.5 if 
a spectral break is present just below the lowest X-ray energy 
detected. GRB afterglows with fi x < fix - 0.5 are classified as 
dark in this picture. 

The complex nature of the X-ray lightcurve revealed by Swift 
between a few tens of seconds and several hours after the trigger, 
places some constraints on the dark GRB identification methods 
mentioned above based on the assumption that both the optical 
and X-ray emission have the same origin, that is, lie on the same 
energy spectrum. The steep initial decay observed in most X-ray 
afterglow light curves and not at optical wavelengths, is indeed 
thought to be produced by the prompt emission, thus by a mech- 
anism different from the one responsible for the afterglow. The 
optical and X-ray spectral energy distributions should then be 
compared after the initial steep decay in order to be sure that the 
prompt emission is not dominating the X-ray afterglow emission 
and that early optical flashes from reverse shocks are not present 
(see, e.g., van der Horst et al. 2009). In addition, in about 70% of 
Swift/XRT detected afterglows, the X-ray light curve between a 
few hundreds of seconds (thus at the end of the steep decay) up to 
several minutes/hours shows a plateau phase that is on average 
not present in the optical counterpart (e.g. Nousek et al. 2006, 
Liang et al. 2009). This behavior is not predicted by the standard 
fireball model. After the plateau phase, alternatively called the 
"shallow phase", the X-ray light curve decays following the af- 
terglow behavior expected from the fireball model. Interestingly, 
the transition from the "shallow phase" to the "normal phase" 
is not accompanied by any spectral variation (see e.g. Liang et 
al. 2009). The interpretation of the so-called "shallow phase" 
observed in the Swift GRB X-ray afterglow remains unclear. 
Therefore, the spectral extrapolation to the optical observation 
times should be interpreted with caution. 



Here we study the 'darkness' properties of two GRBs 
(namely, GRB 100614A and GRB 100615A), which are very 
bright in X-rays, but are not detected in the optical/near-IR band 
and have no reported host galaxy candidate. The paper is orga- 
nized as follows. Sect. 2. summarizes the discovery and observa- 
tions of these two sources, Sect. 3 presents instead our analysis 
method, Sect. 4 illustrates our results, and in Sect. 5 we discuss 
our findings and draw our conclusions. In the following, we as- 
sume a concordance cosmology with Ho = 70 km s _1 Mpc -1 , 
Q m = 0.3, £l A = 0.7. Decay, photon, and spectral indices are in- 
dicated with a, T, and fi, following the standard convention t~ a , 
N~£, and v - ^, respectively. 

2. The dataset 

2.1. GRB100614A 

GRB 100614A was discovered by Swift on June 14, 2010, at 
21:38:26 UT (Stratta et al. 2010). The Swift/BAT (15-350 keV) 
coded mask- weighted light curve shows a relatively smooth peak 
starting at ~ T - 10 s, peaking around T + 50 s, and ending at 
~ T + 275 s, where T is the trigger time. The burst duration has 
been estimated as T90 = 225 ± 55 s, thus classifying this burst 
among the long GRBs. The time-averaged spectrum is best fit by 
a simple power-law model with photon index 1.88 + 0.15. The 
15-150 keV fluence is (2.7 + 0.2) x 10~ 6 erg crrT 2 (Sakamoto et 
al. 2010). 

XRT and UVOT observations were initiated 133 s and 873 
s after the BAT trigger, respectively. XRT immediately found a 
bright afterglow. The astrometrically corrected X-ray afterglow 
position (using the XRT-UVOT alignment and matching UVOT 
field sources to the USNO-B1 catalogue) with 3.4 ks of ex- 
posure, is RA(J2000)=17h 33m 59.82s and Dec(J2000)=+49d 
14' 03.6" (1.7" radius, 90% confidence, Osborne et al. 2010), 
while UVOT did not detect the optical counterpart despite the 
low Galactic reddening in the direction of this burst (E(B-V) of 
0.03 mag, Schlegel et al. 1998). The XRT follow-up (Margutti et 
al. 2010) continued up to about eight days after the trigger and 
the overall light curve shows the canonical "steep-flat-normal" 
decay (e.g. Nousek et al. 2006). 

XMM-Newton also observed GRB 100614A, starting at 
04: 12 UT, on June 15, 2010, for an exposure of 42 ksec (Schartel 
2010), but no results from data analysis have yet been pub- 
lished. Several optical follow-up observations have been per- 
formed with no optical/NIR counterpart detection. Table 1 sum- 
marizes the optical/NIR upper limits to the GRB 100614A af- 
terglow, together with its X-ray flux at 1.7 keV at the corre- 
sponding acquisition times. Magnitudes are in the Vega system 
(unless otherwise stated, see footnotes), but not corrected for 
Galactic extinction. The mean observational epochs from trigger 
are quoted in the second column. Bold face characters indicate 
the epochs with the reddest and deepest upper limits. Fig. 1 (top 
panel) shows the X-ray light curve and all the optical/NIR upper 
limits reported in Table 1 . 

2.2. GRB 1 0061 5A 

GRB 100615A was discovered by Swift on June 15, 2010, at 
01:59:03 UT (D'Elia et al. 2010). The Swift/BAT (15-350 keV) 
coded mask-weighted light curve shows three slightly overlap- 
ping FRED (fast rise exponential decay) peaks. The burst dura- 
tion has been estimated as Tgo = 39 + 2 s. The time-averaged 
spectrum is best fit by a simple power-law model with photon 
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mag 
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mag 


mag 


GRB 1006 14 A" 


6.9 


29.6 


- 


- 


> 18 


- 


- 


- 


- 


- 


GRB 100614A 6 


29 


1.26 


- 


- 


> 21 


- 


- 


- 


- 


- 


GRB 1006 14 A 1 ' 


31 


1.08 


- 


- 


- 


> 21 


- 


- 


- 


- 


GRB 100614A rf 


80 


1.20 


- 


- 


> 24.0 


> 22.8 


- 


- 


- 


- 


GRB 100614A 8 


60 


0.25 


- 


- 


- 


> 24 


- 


- 


- 


- 


GRB 100614A / 


258 


0.10 


>25.2 


>27.1 


>26.4 


>25.9 


>24.9 








GRB 100614A? 


570 


0.07 












>22.7 




>21.6 


GRB 10061 4 A' ! 


1330 


0.05 






> 21.9 












GRB 100615A •• AB 


24 


9.84 




>24.2 


>23.9 


>22.9 


>22.5 


>21.4 


>20.7 


>20.3 


GRB 100615A^ 


28 


9.69 








> 24.0 










GRB 1006 15 A* 


114 


4.92 












> 18.8 


> 17.9 


> 16.6 


GRB 100615A' 


330 


1.97 
















>20.9 


GRB 100615A'" 


13200 


0.013 
















>22.2 



BOOTES-2 TELMA 0.6m robotic telescope data from Jelinek et al. (2010); *: 2-m Liverpool Telescope data from Mundell et al. (2010); c : 2-m 
Liverpool Telescope data from Mundell et al. (2010); d : NOT data from Malesani et al. (2010) e : 4.2m William Herschel Telescope data from 
Levan et al. (2010) f : GTC 10.4-m data from Guziy et al. (2010); s : Gemini-North data from Cenko et al. 2010b; h : AZT-11 telescope data from 
Shakhovskoy et al. (2010); *; GROND data from Nicuesa et al. (2010); j : NTT/Ultracam data from Dhillon et al. (2010); *: PAIRITEL data from 
Morgan et al. (2010); ': Gemini-North data from Cenko et al. (2010a); "': Gemini-North data from Perley et al. (2010); 
AB : Magnitudes given in the AB system 



index of 1.87 + 0.04. The fluence in the 15-150 keV band is 
(5.0 + 0.1) x 10~ 6 ergcirT 2 . 

XRT and UVOT observations were initiated about a minute 
after the BAT trigger. XRT immediately found a bright afterglow. 
The astrometrically corrected X-ray afterglow position (using 
the XRT-UVOT alignment and matching UVOT field sources to 
the USNO-B1 catalogue) with 5.2 ks of exposure is RA(J2000) 
= llh 48m 49.26s and Dec(J2000) = -19d28' 52.4" (1.4" error 
radius at the 90% confidence, Osborne et al. 2010), while UVOT 
did not detect the optical counterpart despite the low Galactic 
reddening in the direction of this burst (E(B-V) of 0.05 mag, 
Schlegel et al. 1998). XRT follow-up continued up to about two 
days after the trigger. As for GRB 100614A, the overall light 
curve shows the canonical "steep-flat-normal" decay. 

A DDT was issued and approved to observe the X-ray af- 
terglow with Chandra (Butler et al. 2010). The source was still 
visible six days after the burst, and the position was refined to 
RA, Dec (J2000) = l\ h 48 m 49 ! .34, -19d 28' 52.0" with uncer- 
tainty of 0.6", and the spectral parameters were found to be in 
agreement with those obtained analyzing the XRT data (see next 
section). The optical/NIR afterglow was searched starting from 
a few minutes from the BAT trigger, using four ground-based fa- 
cilities, but only upper limits to its emission could be set. Table 1 
summarizes the optical/NIR upper limits for the GRB 1006 15 A 
afterglow, together with its X-ray flux at 1.7 keV at the cor- 
responding acquisition times. Magnitudes are in Vega system 
(unless otherwise stated, see footnotes), but not corrected for 
Galactic extinction. The mean observational epochs from trigger 
are quoted in the second column. Bold face characters indicate 
the epochs with the reddest and deepest upper limits. Fig. 1 (bot- 
tom panel) shows the X-ray light curve and all the optical/NIR 
upper limits for GRB 1 006 1 5 A. 



3. Data analysis 

We extract the optical to X-ray spectral energy distribution 
(SED) for both GRBs by selecting those epochs at which we 
have the deepest and reddest observations (i.e. less affected by 
any dust extinction), to constrain at best the intrinsic optical af- 



terglow flux upper limit. In addition, we attempt to select an 
epoch not too close to the initial X-ray steep decay, which is 
thought to be produced by a different component than the one re- 
sponsible for the afterglow emission. Magnitudes were corrected 
for Galactic absorption. 

Swift/XRT data were calibrated, filtered and screened us- 
ing the XRTDAS package included in the HEAsoft distribu- 
tion (v6.10) as described in the XRT Software User's Guid^B 
Unabsorbed X-ray fluxes were estimated at the selected epochs 
(see below). 

A broken power-law model was fitted to the data following 
the van der Horst et al. (2009) method, therefore fixing the SED 
normalization and the high-energy spectral index to the value ob- 
tained from our X-ray data analysis (within its 90% confidence 
range), the spectral break at the X-ray energies and the optical to 
X-ray spectral index as Pox — fix - 0.5. 

We model the optical suppression from the X-ray extrapola- 
tion assuming either a Milky Way (MW) or a Small Magellanic 
Cloud (SMC) extinction curve. We also test the attenuation curve 
obtained for a sample of starburst galaxies (Calzetti et al. 1994). 
We consider the upper limits as positive detections, and we ver- 
ify that the model-predicted fluxes are consistent with the data, 
i.e. equal to or below the upper limits. 

3.1. GRB100614A 

For this burst, the X-ray steep decay phase ends about 30 min- 
utes after the trigger, while the plateau stops 2.27 days after 
the trigger. The reddest and deepest upper limits are those ob- 
tained 258 and 570 minutes after the trigger with the 10.4-m 
Gran Telescopio Canarias (GTC) in the ugriz photometric bands 
(Guizy et al. 2010) and with Gemini Near InfraRed Imager on 
the 8-m Gemini North telescope in the J and K bands (Cenko et 
al. 2010b), respectively (see Table 1). The data of the first ob- 
servation were photometrically calibrated using the SDSS stars 
and the upper limits are at the 3cr confidence level (Guziy et al. 
2010; Guziy, private communication). The data of the second ob- 
servation were calibrated using three 2MASS stars and are given 

1 http://swift.gsfc.nasa.gov/docs/swift/analysis/ 
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Fig.l. The GRB 1006 14A (upper panel) and GRB 100615A 
(bottom panel) X-ray lightcurves, together with all the opti- 
cal/NIR upper limits estimated using several ground-based fa- 
cilities (see Table 1 for details). Cyan, blue, magenta, green, yel- 
low, black and red upper limits represents the g' (and all those 
bands bluer than g'), R, I, z, J, H, K bands, respectively. Vertical 
dashed lines emphasize the reddest observations and the time at 
which they have been obtained. Such observations are used in 
our analysis. The faintest X-ray data in the bottom panel draws 
the GRB 100615A Chandra observation (Butler et al. 2010). 
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Fig. 2. GRB 100614A SEDs at two different epochs: GTC ugriz 
upper limits taken 4.3 hours after the trigger (upper panels) and 
the NIR J and K band taken 9.5 hours after the trigger (lower 
panels), with fiox fixed at fix - 0.5 (solid red line) and the op- 
tical suppression modelled with the MW or starburst extinction 
curves. Black dotted line represents the best-fit value for the X- 
ray spectral slope (fix = 1 -50). Dashed red lines enclose the 90% 
uncertainty in fiox- 



ing to the magnitudes of the first observation are the following 
(from the U to the z band): f v < 0.13, f g < 0.056, f r < 0.085, 
fi < 0.1 1, f z < 0.24 yuJy. The optical upper limits corresponding 
to the J and K data points of the second observation are fj< 1.3 
liiy and f K < 1 .5 yuJy, respectively. 

The 0.3-10 keV energy spectrum integrated between 10 ks 
and 100 ks after the trigger, that is during the plateau phase, 
is closely fitted by an absorbed power law with an equivalent 
hydrogen column density Nh — (1.6 ± 0.7) x 10 21 cirT 2 beyond 
the Galactic one measured at the position of the X-ray afterglow 
(N H ,cai = 2.2 x 10 2() cm" 2 , Kalberla et al. 2005) and a photon 
index of F = 2.5 + 0.3 (90% confidence range). The count rate 
to unabsorbed flux conversion factor is 6 x 10 _1 1 erg crrT 2 cts -1 . 

The X-ray spectrum does not vary significantly during the 
plateau phase, thus we assume this spectral shape for the epochs 
at which we wish to extract the broad-band SEDs, and we com- 
pute the X-ray normalization level from the light curve at the 
given epochs. The count rate 258 and 570 minutes after the trig- 
ger is c\ = 0.03 counts s _1 and c 2 = 0.02 counts s _1 . Converting 
these values into unabsorbed flux using the conversion estimated 
factor, we find /i = 1 .8 x 10~ 12 erg cirT 2 s _1 and f 2 = 1 .2 x 10~ 12 
erg cirT 2 s _1 . The flux densities at 1.7 keV (logarithmic mean of 
the XRT energy range) are therefore f Vl = 0.10 and f Vl = 0.07 
//Jy 258 and 579 minutes after the trigger, respectively. 



at the 5cr confidence level. The RMS spread evaluated by com- 
paring these data with calibrators is 0. 1 (0.2) mag in the J (K) 
band (Cenko et al. 2010b; Cenko, private communication). All 
the available magnitudes were corrected for Galactic reddening 
in the direction of this burst. The optical upper limits correspond- 



3.2. GRB 10061 5 A 

For this burst, the X-ray steep phase ends 185 s after the trigger, 
and the flat one about 5 ks later. The reddest observations for the 
GRB 1006 15 A afterglow were performed in the K band 20 min- 
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utes (Nicuesa et al. 2010), 5.5 hours (Cenko et al. 2010a), and 
9.2 days (Perley et al. 2010) after the BAT trigger (see Table 1). 
The first observation falls during the plateau phase, while the last 
two take place during the normal phase. The magnitudes of the 
first observation are obtained using the GROND zero point and 
the 2MASS catalog, and are at the 3cr confidence level (Nicuesa 
et al. 2010). The data of the second observation are calibrated 
against three 2MASS stars and the upper limits are at the 5<x 
level. The RMS spread when comparing these data with cali- 
brators is 0.2 in the K band (Cenko et al. 2010a, Cenko, private 
communication). Finally, the data of the last observation are cal- 
ibrated using two 2MASS stars and the upper limits are at the 3<x 
level. The calibration is accurate to 0.20 - 0.25 mag (Perley et 
al. 2010, Perley, private communication). All these magnitudes 
were corrected for Galactic reddening in the direction of this 
burst. The optical upper limits corresponding to the magnitudes 
of the first observation are the following (from the g' to the K 
band): f g , < 1.0,/> < 1.0, /? < 2.3, f t < 2.8, fj < 9.56, 
fn < 18.7, fx < 27.9 /vJy. The optical upper limits correspond- 
ing to the K data points of the second and third observations are 
fx < 16.1 yt/Jy and fx < 4.9 /iJy, respectively. 

The 0.3-10 keV energy spectrum was integrated during the 
plateau, the normal and the normal+plateau phases. In all cases, 
it could be closely fitted with an absorbed power law model. The 
photon indices and Ng values are consistent in the three integra- 
tion intervals, thus we assume as the most accurate estimates of 
the spectral parameters, those for the overall integration, which 
are F = 2.35 + 0.15 and N H = (1.05 + 0.12) x 10 22 cirT 2 be- 
yond the Galactic value at the position of the X-ray afterglow 
(NH Ca , = 3.3 x 10 20 cnr 2 , Kalberla et al. 2005), with a count 
rate to unabsorbed flux conversion factor of 1 .5 x 10~ 10 erg cm -2 
cts~'. 

As for GRB 100614A, we assume this spectral shape for the 
broad-band SED analysis, computing the X-ray normalization 
level from the light curve at the give epochs. The count rate at the 
selected epochs is c\ = 1.0 cts/s, ci = 0.2 counts s _1 , and cj, = 
0.005 counts s . Converting these values into unabsorbed flux 
using the conversion estimated factor, we find f\ = 1.5 x 10~ 10 
,f 2 = 3.0x10"' 1 erg cnr 2 s -1 , and f 3 = 7.5 x 10~ 13 



erg cnr 2 s _1 , 

erg cnr 2 s _1 . The flux densities at 1.7 keV (logarithmic mean of 
the XRT energy range) are therefore / V] = 9.84, / y , = 1.97, and 
./:,. -0.013//.I.V. 



4. Results 
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Fig. 3. GRB 1006 15A SEDs at three different epochs: GROND 
grizJHK upper limits taken 24 minutes after the trigger (upper 
panels) and Gemini-North K and Ks upper limits taken 5.5 hours 
and 9.2 days after the trigger, respectively (middle and lower 
panels), with Pox fixed at fix - 0.5 and the optical suppression 
modelled with the MW or starburst extinction curves. Black dot- 
ted line represents the best-fit value for the X-ray spectral slope 
(fix = 1.35). Dashed red lines enclose the 90% uncertainty in 
Pox- 



We find that GRB fluxes computed from both optical and NIR 
data are well below the most conservative extrapolation from X- 
rays and require strong absorption using the data taken within 1 
day after the trigger. For GRB 1006 15 A, very late-time data are 
available (9.2 days after the trigger), for which the NIR flux is 
still below the X-ray extrapolation. 

The GRB 100614A SED with the reddest flux upper limit 
(i.e. 570 minutes after the trigger) requires a rest-frame V-band 
dust extinction of Ay > 47, 58, and 11 mag, assuming a MW, 
SMC, or a starburst extinction curve, respectively, and using the 
best-fit value for fix (Fig. 2). Even fixing the optical to X-ray 
energy spectral index to its lowest allowed value (within its 90% 
confidence range), that is, in the most conservative case, results 
still provide very high Ay lower limits (Table 2). Less stringent 
constraints on Ay are obtained using the GTC ugriz flux upper 
limits obtained 258 minutes after the trigger. We obtain Ay > 13 
mag with either the MW and SMC extinction curve and Ay > 8 
mag with the starburst case. 



For GRB 100615A, we obtain even tighter lower limits. The 
SED extracted 24 minutes after the burst requires Ay > 64, 79, 
or 16 mag assuming either a MW or SMC extinction curve, or 
a starburst attenuation curve, respectively. These lower limits 
are still very high for the SED extracted 5.5 hours post burst: 
Ay > 58, 72, and 15 mag for the three extinction recipes. Less 
critical but still high values, are obtained even 9.2 days from the 
burst: Ay > 22, 27, and 6 mag (Fig. 3). The latter lower limits are 
lower than the ones obtained at earlier epochs (but still extreme), 
possibly due to a selection effect, since at later times the X-ray 
flux decreases, but the optical/NIR upper limits can not become 
fainter consistently, owing to the instrument detection limits. 

Table 2 reports all the Ay lower limits evaluated from the 
GRB 100614A and GRB 100615A data, for MW, SMC, and 
starburst extinction laws at the given mean observation epochs. 
These lower limits are computed for the reported ranges of the 
X-ray spectral index (fix = T — 1, estimated in Sect. 3) that cor- 
responds to the minimum, maximum, and mean value of the es- 
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Table 2. Lower limits to the visual extinction towards 
GRB 100614A and GRB 100615A. 



Source 


t-T 


fix 


Pox 


Av,MW 


Av,SMC 


Av.sb 




min 






mag 


mag 


mag 


GRB 100614A 


258 


1.20 




8 


8 


6 


GRB 100614A 


258 


1.50 


< 0.1 


13 


13 


8 


GRB 100614A 


258 


1.80 




18 


18 


11 


GRB 100614A 


570 


1.20 




25 


31 


6 


GRB 100614A 


570 


1.50 


< 0.4 


47 


58 


11 


GRB 100614A 


570 


1.80 




69 


85 


17 


GRB 1006 15 A 


24 


1.20 




39 


47 


10 


GRB 100615A 


24 


1.35 


< 0.1 


49 


60 


13 


GRB 1006 15 A 


24 


1.50 




59 


73 


16 


GRB 1006 15 A 


330 


1.20 




47 


58 


12 


GRB 100615A 


330 


1.35 


< 0.0 


58 


72 


15 


GRB 1006 15 A 


330 


1.50 




69 


86 


18 


GRB 1006 15 A 


13200 


1.20 




12 


15 


3 


GRB 100615A 


13200 


1.35 


< 0.5 


22 


27 


6 


GRB 1006 15 A 


13200 


1.50 




33 


40 
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timated 90% confidence range. Bold face characters indicate the 
results obtained with the most probable X-ray spectral index. 
Upper limits to the optical-to-X-ray spectral indices Pox are also 
shown. 



5. Discussion 

We have analyzed two Swift 'dark' GRBs, namely, 
GRB100614A, and GRB 1006 15 A. These GRBs are dark 
according to every definition proposed until now. They are not 
detected in the optical/NIR down to very faint limits, despite 
follow-up campaigns at ground-based facilities began a few 
minutes from the BAT triggers (see Table 1). In addition, their 
optical-to-X-ray spectral indices satisfy Pox < 0.5 (Jakobsson et 
al. 2004 criterion) and /3 x < fix - 0.5 (van der Horst et al. 2009 
criterion). The identification of these two GRBs as dark bursts 
according to the above-mentioned methods, is the consequence 
of their intense X-ray flux coupled to the optical/NIR missing 
detections. GRB 100614A (GRB 100615A) indeed falls in the 
upper 30% (5%) of the distribution of SwiftPiRT GRB fluxes at 
1 1 hr from the burst (Gehrels et al. 2008). 

The outcome of our analysis is surprising. To explain the 
deepest NIR upper limits (i.e. the less affected by dust) in 
terms of a flux suppression described by either a MW or SMC 
dust extinction laws, A v > 47 (A v > 58) mag is needed for 
GRB 100614A (GRB 100615A) before one day and A v > 22 at 
9 days after the trigger for GRB 1006 15 A. Such extreme Ay val- 
ues have never been observed before and require an explanation. 

As an example, Perley et al. (2009) studied a sample of 29 
Swift bursts rapidly observed by the Palomar 60-inch telescope 
(14 of which were classified as dark) from April 2005 to March 
2008. From both optical to X-ray afterglow spectral analyses 
and host galaxy studies (applying the p x < 0.5 criterion by 
Jakobsonn et al. 2004 and assuming in general a SMC extinc- 
tion law), they find that more than half of the dark sample is 
extincted by dust, with three events featuring A v > 2-6 mag, 
two GRBs possibly being high redshift events (4.5 < z < 7), 
and three being underluminous events (not detected in the op- 
tical bands, but not dark according to Jakobsonn et al. (2004) 
criterion). The dust extinction in their sample is compatible to 
that observed along the dustiest MW sightlines. Extinction val- 



ues of Ay < 5 are found along more than 700 Galactic sight- 
lines (Diplas & Savage 1994). In a handful of extreme cases, Ay 
can be higher, up to ~ 20 (Predehl & Schmitt 1995). Rol et al. 
(2007) studied what they called a 'prototype' dark GRB, namely, 
GRB 051022. This burst misses any optical/NIR detection, such 
as GRB 100614A and GRB 100615A, and is the darkest burst 
in the paper by van der Horst et al. (2009) according to their 
classification criterion. Applying our analysis to GRB 051022, 
we find Ay > 16 and 20 mag assuming a MW and SMC ex- 
tinction curve, respectively. Again, these values are well below 
the Ay lower limits we need to fit the deepest NIR upper limits 
of GRB 100615A and GRB 100614A, using the same extinction 
curves. 

While a SMC-like extinction curve can adequately fit a large 
fraction of the dust extinction from GRB host ISM, our present 
picture of GRB host galaxies makes an extremely obscured en- 
vironment of this kind a very unlikely possibility. Indeed, GRBs, 
even reddened ones, are hosted by blue or normal galaxies, 
which are commonly detected in ordinary galaxy surveys. This 
favours a scenario of a host morphology where the line of sight 
to the GRB is dusty, i.e., dust obscures only localized regions 
(see e.g., Perley et al. 2009 and references therein). An in situ 
obscuration appears to be insufficient to be responsible for the 
extreme extinction levels we measure. 

All these considerations hold if the optical radiation and X- 
rays are part of the same synchrotron spectrum. They could orig- 
inate from different emission processes or even be produced in 
different, independent emission regions. This is possible in par- 
ticular during the so-called "shallow phase", where X-ray emis- 
sion may be dominated by an emission component that differs 
from the one from which the optical flux originates (e.g. Zhang 
et al. 2006 for a review). Interesting hints about the early optical 
to X-ray afterglow spectral behavior came from Greiner et al. 
(201 1). These authors addressed the darkness problem by study- 
ing all the GRBs observed by the GROND imager mounted on 
the 2.2 m MPI/ESO telescope at La Silla (Greiner et al. 2008). 
Broad-band SEDs of those GRBs with simultaneous optical and 
X-ray data within the first 240 minutes and after the early X-ray 
steep decay (43 GRBs, 39 of which are long) have been extracted 
and a simple or a broken power-law spectral model has been fit- 
ted to the data according to the darkness criterion of van der 
Horst et al. (2009). The optical to X-ray SEDs of those GRBs at 
known redshifts clearly show that in no case the optical fluxes 
are above the X-ray extrapolation, a signature that would have 
confirmed the possible distinct origin of the early X-ray emis- 
sion from the optical one. In particular, modelling the optical 
suppression from the X-ray extrapolation as SMC- or MW-type 
dust extinction, Greiner et al. (2011) found in all cases a good 
agreement between the expected fireball spectral model and the 
data, and a fraction of 25% - 40% of the bursts of their sample 
were found to be dark (different percentages depend on the defi- 
nition used), where their darkness can be explained by moderate 
extinction (0.5 < A v < 1.5), or high redshift (z > 5) for 22% 
of the dark bursts. These results show the lack of an evident in- 
consistency between the optical and X-ray early afterglow SED. 
However, the brightness of the X-ray flux may be biased by the 
possible presence of a dominant component that differs from the 
one responsible for the optical emission: excess of the X-ray to 
optical flux ratio could thus mimic a stronger optical absorp- 
tion than the real one. Since the two SEDs of GRB 100614A, 
the first one of GRB 1006 15 A, and (marginally) the second one 
of GRB 1006 15 A are all extracted during the shallow phases of 
these GRBs, a different origin of the optical and X-ray emission 
in these epochs could at least in part explain the optical darkness 
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of our GRBs. However, we note that for GRB 100615A we have 
extracted a SED at a very late time, about nine days after the 
end of the plateau phase, and we have still obtained very high 
Ay lower limits (A v > 20) assuming either a MW or a SMC ex- 
tinction curve. These values are less extreme than that obtained 
using the other SEDs, but still very high, suggesting another or at 
least a concurring mechanism to account for the GRB darkness. 

This complementary explanation could be that local extinc- 
tion recipes, such as MW or SMC ones are inadequate for re- 
producing the optical suppression in the host galaxies of these 
two GRBs. Indeed, the dust and extinction properties of GRB 
host galaxies are still poorly known. For example, modelling 
the dust absorption using greyer extinction laws, such as the 
attenuation curve obtained from the observations of starburst 
galaxies proposed by Calzetti (1994), brings GRB 100614A and 
GRB 100615A to require Ay lower limits that are less extreme, 
despite still being very high. More gray extinction curves have 
already been invoked for GRB environments by several authors 
(e.g. Perley et al. 2009; Stratta et al. 2004, 2005; Chen et al. 
2006; Li et al. 2008). The shape of the extinction curve provides 
information on the dust properties. The study of dust in GRB 
environments is particularly useful for studying any evolution of 
dust properties and dust production mechanisms up to very high 
redshifts and in "normal" galaxies. The latters are indeed more 
representative of the majority of the galaxy population, rather 
than the more extreme galaxies, such as those hosting quasars. 

A mixture of moderate-to-high redshift and extinction can 
reduce the dust level necessary to explain the SEDs. In fig. 4, we 
plot as an example the Ay values as a function of z obtained for 
the second epoch of GRB 100615A observations. Although the 
visual extinction is considerably lower than for the z — case, 
Ay ~ 10 at z — 2 and Ay ~4-5atz = 5is still required, 
regardless of the adopted extinction recipe. Kann et al. (2010) 
reported the computed rest-frame Ay values for a vast sample of 
GRBs with known redshift (see their fig. 3). All GRBs in that 
sample are modelled with Ay < 1.3. Similar results are obtained 
by Greiner et al. (2011) whose SEDs of GBRs at known red- 
shift are modelled with Ay < 1.5. We stress that GRB samples 
with known redshift are usually biased towards non-absorbed 
GRBs, and the comparison between their Ay and that of our 
sources must then be taken with caution. However, the Greiner 
et al (201 1) sample is more than 90% complete, since only 3 out 
of the 39 long GRBs observed by GROND miss a redshift esti- 
mation. Although our Ay lower limits in this scenario are well 
above those estimated for the Kann et al. (2010) and Greiner 
et al. (2011) samples, they are not extremely unlikely. Thus, an 
intrinsic origin and/or dust extinction, coupled to a moderately 
high redshift could explain the darkness of GRB 100614A and 
15A. 

A more exotic but intriguing possibility would be that these 
GRBs are extremely high redshift events. Assuming that the lack 
of any detection in the reddest NIR band (K-band) is due to Lya 
absorption from the intergalactic hydrogen neutral fraction, we 
can set a redshift lower limit of z > 17. The first population of 
very massive stars (PopIII, 1O 2 M < M < 1O 3 M ) is expected 
to form at z ~ 20. Their death is supposed to leave behind black 
holes of several tens of solar masses, which could be the early 
progenitors of active galactic nuclei. These fast-spinning black 
holes have a rotational energy of ~ 10 55 erg or more that can 
power a GRB explosion (see e.g., Komissarov & Barkov 2010; 
Meszaros & Rees 2010; Suwa & Yoka 201 1). The isotropic ener- 
gies of GRB 100614A and GRB 100615Acanbe estimated from 
their BAT fluence (Sakamoto et al. 2010; Palmer et al. 2010). 
The resulting values assuming z - 18 are 1.3 x 10 54 erg and 




Fig. 4. The visual extinction Ay as a function of redshift for the 
second epoch SED of GRB 1006 15 A. Green solid, red dotted, 
and blue dashed lines are for the Milky Way, Small Magellanic 
Cloud, and starburst extinction recipes, respectively. 



7.2 x 10 53 erg, respectively. Such energies fall in the bright tail of 
the GRB distribution, the most energetic burst detected to date, 
GRB080916C, having E iso = 9 x 10 54 erg (Abdo et al. 2009). 
The Amati relation (Amati et al. 2008) allows for a 2<x scat- 
tering of the prompt emission energy peak of 0.5 - 5 MeV in 
the rest frame, assuming an isotropic energy of ~ 10 54 erg. The 
reported isotropic energies must be considered as lower limits 
both because of the conservative choice of the z used and be- 
cause the BAT detector does not constrain the position of the 
peak emission preventing a bolometric estimate of the emitted 
energy. Thus, that these peaks are not required in the 15-150 
keV spectral fit of the GRB 100614A and GRB 100615A BAT 
data (Sakamoto et al. 2010; Palmer et al. 2010) is not surprising. 
In addition, according to some models (see e.g., Suwa & Yoka 
201 1), a considerable fraction of the available energy to produce 
the GRB must be used to pierce the envelope of the very massive 
star. The lack of detection of any host galaxy candidate for these 
GRBs represents additional support of this scenario. On the other 
hand, against the high redshift interpretation there is that the ex- 
pected GRB rate at z > 17 is extremely low, between 0.5 and 1 
GRB every 10 yr (Bromm & Loeb 2006). In addition, there is the 
non-negligible equivalent hydrogen column density value mea- 
sured for GRB 100614A and GRB 100615A from X-ray spec- 
troscopy. For example, at z = 18, the rest-frame column density 
would be roughly ( 1 +z) 8 ^ 3 times higher than the measured value, 
thus ~ 300 - 3000 x 10 22 cnr 2 for our GRBs. These high hydro- 
gen column densities would require an extremely dense environ- 
ment surrounding the GRB, which challenges this scenario. 

A possible way to differentiate between extremely high- 
redshift, exotic extinction recipes and emission from distinct 
components for sources such as these, would be to search for 
the afterglow in the mid- or far-IR bands. A non-detection also 
in these bands could hardly be explained using any extinction 
law and would definitely rule out a high redshift origin. A multi- 



8 



V. D'Elia and G. Stratta: GRB 100614A and 100615A: two peculiar dark GRBs 



band detection compatible with the X-ray flux at late times (e.g. 
> 10 5 s after the trigger, i.e. after the end of the "shallow phase") 
would instead favour the fireball model. In this case, given the 
lack of NIR detections, the origin of the darkness would be in ex- 
otic extinction or high redshift, depending on the spectral shape 
in the mid- and far-IR bands. 
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